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Abstract. We numerically investigate under which conditions the planet detected at 2.1 AU of 7 Cephei could form 
through the core-accretion scenario despite the perturbing presence of the highly eccentric companion star. We first 
show that the initial stage of runaway accretion of kilometer-sized planetesimals is possible within 2.5 AU from the 
central star only if large amounts of gas are present. In this case, gaseous friction induces periastron alignment of 
the orbits which reduces the otherwise high mutual impact velocities due to the companion's secular perturbations. 
Fhe following stage of mutual accretion of large embryos is also modeled. According to our simulations, the giant 
impacts among the embryos always lead to a core of 10 M$ within 10 Myr, the average lifetime of gaseous discs. 
However, the core always ends up within 1.5 AU from the central star. Either the core grows more quickly in the 
inner region of the disc, or it migrates inside by scattering the residual embryos. 
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1. Introduction 

Among the 15 presently known binary star systems har- 
bouring extra solar planets, 7 Cephei is that with the 
closest companion star (with the exception of Gliese 86, 
where the possible co mpanion is be lieved to be a brow n 
dwarf lEls et all 12001^ . According to lHatzes et all l|2003|) . 
the secondary star has an orbit with a semimajor axis of 
18.5 ± 1.1 AU and an eccentricity of 0.361 ± 0.023. The 
planet detected around the primary K0 III giant star, has 
a mass of Msin i = 1.7±0.4 Jupiter masses and an orbital 
semimajor axis of 2.13 AU. The planet is located inside a 
stable region with a dynamical lifetime of at least 1 Gyr. 
This region extends to about 4 AU from the cen tral star 
l|Holman and Wiegertlll999HDvorak et alll2003|h 

Giant planets orbiting one of the stars in a binary sys- 
tem offer the possibility to test the core-accret ion model 
for giant planet formation l|Pollack et all Il996[) in a com- 
plex dynamical environment. The vicinity of a companion 
star on a highly eccentric orbit may prevent planetary for- 
mation because the companion reduces the size of the ac- 
cretion disc, and it excites high relative velocities between 
colliding planetesimals. 

In this paper we study the in situ formation of the 
giant planet in the 7 Cephei system by numerically in- 
vestigating both the stage of planetary embryo accretion 
from planetesimals and the following stage of high veloc- 
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ity collisions between large embryos which form the core 
of the giant planet. This last stage precedes that of rapid 
gas accretion which sets on after sufficient accumulation 
of mass onto the core. 

A critical parameter for the initial stage of planetesi- 
mal accretion is the relative impact velocity Av. This ve- 
locity determines whether accretion or erosion dominates 
the planetesimal collisional evolution. In a binary star sys- 
tem, the secular perturbations of the closeby companion 
star may play a critical role by exciting the relative veloc- 
ities among the planetesimals. Too large Av prevent plan- 
etesimals from growing by accretion and the initial plan- 
etesimal population may even be ground down to dust. We 
explored the distribution of collisional relative velocities of 
a planetesimal population surrounding the primary star of 
the 7 Cephei system by using a deterministic code that in- 
cludes the effect of gas drag. This later mechanism might 
indeed play a crucial role since gas drag might force peri- 
astron and eccentricity alignment of the perturbed plan- 
etesimals orbits and thus reduce the large Av induced by 
the c ompanion's perturbations (e.g. iMarzari and Scholl . 
120001) . 

We also investigated the late stage of core formation 
by following a population of large planetary embryos with 
a numerical model ba sed on Chambers' Mercury code 
llChambers et all l2002j) which takes into account all mu- 
tual gravitational forces among the embryos as well as mu- 
tual collisions. Our main concern is here to see if mutual 
embryos accretion can lead to a final planet at the right 
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Fig. 1. Distribution of particles eccentricity vs. semimajor 
axis after t = 50000 years for a simulation where the gas 
drag is not included. 

place, i.e. ~ 2.1 AU, and within 10 7 year s, i.e. the typical 
survi val time for circumstellar discs fe.g. lArmitage et all 
120031) . 

The dynamical conditions for accretion within a plan- 
etesimal swarm perturbed by the companion of 7 Cephei 
and affected by gas drag are investigated in section 2. 
Section 3 is devoted to the study of the final accumulation 
of massive embryos into a core. In Section 4 we discuss the 
results. 

2. Relative velocities among planetesimals 

Planetesimal accretion in a circumstellar disc occurs when 
mutual encounter velocities are lower than the escape ve- 
locity Vesc of the mutually colliding bodies (corrected by 
a factor that accounts for the energy dissipation in the 
impact). In the absence of external perturbations, average 
eccentricities and inclinations in a planetesimal swarm are 
small, and relative velocities are always low enough to al- 
low fast runaway accretion of lar ger bodies at tim escales 
of the order of 10 4 -10 5 years fe.g iLissauerl Il993j) . When 
the planetesimal population surrounds a star in a binary 
system, the gravitational pull of the companion star ex- 
cites large eccentricities and the impact velocities may in 
some regions exceed by far the mutual escape velocity. In 
these regions the swarm would erode to dust rather than 
form a planet. 

A first step to test whether planetary formation is 
possible around 7 Cephei is thus to estimate the relative 
velocities between planetesimals orbiting within the dy- 
namically stable region around the star. We have used 
a numerical code alre ady adopted in similar studies (see 
iThebault et all |2002l and references therein). This code 
computes the orbits of a swarm of massless particles under 
the influence of one (or several) gravitational perturbers. 
The particles can also be subjected to gas drag forces. 
During the orbital computations, the relative velocities 



Fig. 2. Plot of Av vs. radial distance to the star of the 
planetesimal swarm at t = 50000 years. Each dot stands 
for an encounter that has occurred at a velocity Av and 
at a distance r from the star during [t-10 4 ,t+10 4 ] . The 
continuous line is the average encounter velocity at each 
distance from the star. 



during all the mutual encounters are recorded. These ve- 
locities approximate with a satisfying degree of precision 
the impact velocities during planetesimal collisions. Their 
distribution tells us whether planetesimals accrete into 
larger bodies or erode into smaller pieces. 

In all our simulations we integrate the orbits of 2500 
test particles initially distributed in the 0.3-5 AU region. 
The initial eccentricities and inclinations satisfy the re- 
lation i = e/2 and are chosen such that the average en- 
counter velocity is (Aw) ~ lOm.s^ 1 . We stop our simula- 
tions after tfi na i = 10 years, a typical timesca le for the 
formation of planetary embryos <|Lissaueii ll993'l. 

2.1. Gravitational model without gas drag 

We performed a first test simulation where only the grav- 
itational forces of the two stars are accounted for. Fig^ 
shows the major secular perturbations of the companion 
star on the planetesimal population. At the beginning of 
the simulation eccentricities progressively increase but the 
periastra of these forced orbits are almost aligned. As time 
goes on, however, there is a progressive dephasing pro- 
cess due to the different values of semimajor axis of the 
planetesimals in the swarm. This dephasing leads to the 
wavy pattern observed in Fig^ The planetesimals orbit- 
ing close to the outer limit of the disc at ~ 4.5 AU are more 
strongly perturbed by the secondary star and they reach 
very high eccentricities. They also survive over 10 5 years 
and they cross the orbits of inner planetesimals in the 
3AU region at high Av (hence the wide "wing" of high 
Av that is superimposed on the main swarm of lower rel- 
ative velocities in Fig|2J). The contribution of these scat- 
tered planetesimals is critical since they increase relative 
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Table 1. Initial parameters for the reference gas drag runs 
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Fig. 3. Evolution of the average encounter velocity distri- 
bution, at 5 different epochs, for the gravitational model 
without gas drag in a disc initially truncated at 4AU. 
Race is the size of the smallest body that can accrete 
matter, after a collision with an equally sized object, for 
a given encounter velocity Av. It is given by R aC c = 
( \ ' 5 

/ ( s „p pG J Av, where p p = 3 g.crn" 3 is the density of the 
body and / is a coefficient accounting for the energy dissi- 
pation in an impact. We adopt here the usual assumption 
that 90% of the energy is dissipated in the radial direction 
(e.g. Petit and Farinella 1993) which leads to (/} ~ 0.7 



velocities far beyond the accumulation threshold every- 
where in the inner disc. However, it is a matter of debate 
whether planetesimals could indeed form in the outer edge 
of the stability region around the primary star. Moreover, 
theoretical calculations of binary-disc interactions pre- 
dict that companions might truncate circumstellar discs 
at an o uter radius of 0.2-0.5 times the binary semima- 
jor axes l|Artvmowicz and Lubowiri994|) . The grain coag- 
ulation process and the first impacts between the proto- 
planetesimals would in any case have been highly ener- 
getic, because of the companion's perturbations, possibly 
halting any further growth. For this reason, we performed 
an additional simulation with the planetesimal disc trun- 
cated beyond 4 AU. In Fig|21we show the distribution of 
the impact velocities at different evolutionary times in this 
truncated disc. At the beginning, the secular oscillations 
do not induce a significant increase of Av among plan- 
etesimals because of the strong phasing of all neighbour- 
ing orbits. Nevertheless, the oscillations get narrower with 
time and e and periastron uj gradients between adjacent 
regions keep increasing. At some point, the orbital phas- 
ing is no longer strong enough to prevent orbital crossing 
between bodies with different semimajor axis. The inner 
limit for the region where orbit crossing occurs is very 
sharp, with (Av) increasing from lOm.s" 1 to ~ lkm.s -1 
in less than 0.2 AU (FigO. Furthermore, this limit evolves 
inwards with time, so that the 2.1 AU region is reached in 
less than 1.5 x 10 4 years (FigOJ. As a consequence, the 
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timespan during which accretion of km-sized objects is 
possible within 2.1 AU from the star, assuming that the 
swarm was initially aligned, is typically of the order of 
10 4 years. Numerical accretion models predict that this 
time span is in principle enough to allow the formation 
of 100 to 1000 km-sized bodies (e.s. IWetherill fc Steward. 
ll993UBarge fc Pellatl Il99i IWeidenschilling et all Il997l) . 
However, the large relative velocities that build up after 
10 4 years might halt any further growth of the accreted 
objects. It is even possible that the planetesimal forma- 
tion process from the dust of the circumstellar disc does 
not necessarily lead to orbits that are initially aligned. In 
this case, the planetesimals would have from the begin- 
ning of their evolution high relative velocities. This might 
be the case if we suppose that the binary forms by direct 
stellar-like gravitational instabilities. The companion star 
would reach its present mass well before the onset of plan- 
etesimal accretion in the inner disc. Thus, the formation 
of large embryos in the 2 AU region probably requires the 
presence of some additional mechanism. 

2.2. The effects of gas drag 

Frictional drag by the gas of the protoplanetary disc is an 
important factor in early planetary formation. It affects 
planetesimal orbits in two ways: 

— It restores the periastron alignment 

l|Marzari and Schorl l200fl) . preventing orbital 
crossing of orbits with different semimajor axes. At 
the same time, it partially damps the amplitude of 
oscillations in eccentricity induced by the companion 
star. 

— It causes a drift towards the central star that is par- 
ticularly fast for planetesimals in binary star systems: 
The forced component in the planetesimal eccentricity 
is large, in spite of the damp effect, and it causes a fast 
drift towards the star. 

Following lWeidenschilling and Davisl l)l985|) . we model 
the gas drag force as: 

F = -Kvv (1) 

where F is the force per unit mass, v the velocity of the 
planetesimal with respect to the gas, v the velocity mod- 
ulus, and K is the drag parameter. It is a function of the 
physical parameter of the system and is defined as: 

K = P^- (2) 
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Fig. 4. eccentricities vs. semi-major axis at t = 2500 years 
for 2 00m-sized planetesimals. The gas density is taken 
from iBodenheimer et alJ l(2000h with a nominal Hayashi 
(1981) radial profile. 



where p g is the gas density, p p i and R p i the planetesi- 
mal density and radius, respectively. Cd is a dimension- 
less coefficient related to the shape of the body (~ 0.4 
for spherical bodies). As it appears from the expression 
of the parameter K, the relevance of the drag force on 
the dynamics of the planetesimals depends on the ratio 
between the gas density and the planetesimal size. As a 
consequence, the results we obtain for a particular value 
of p g and R p i can be extended to different combinations 
of the two parameters without any additional simulation. 

We take as a reference va lue for the gas density that 
of IBodenheimer et all l|2000f) who modeled the in-situ 
formation of a giant planet around 47UMa at 2AU. 
Of course, the two systems are not exactly comparable, 
the planet around 47UMa being ~ 50% more massive 
whereas the star is ~ 60% less massive. While these two 
values might compensate one another it is difficult to 
quantitatively estimate ho w much exactly. However, the 
IBodenheimer et all J2000) study gives a good reference 
value for the in-situ formation of a giant planet on an inner 
orbit. According to these authors, this formation requires 
a local gas density p g (2AU) — 2 X 10 _9 g cm -3 , about an or- 
der of magnitu de higher than the value deduced from the 
lHavashil l)l981^ minimum mass solar nebula. We here use 
this value for p g (2AU) an d adopt the classical p g oc r~ 7 
radial profile of Hayashi (1981). Different values of R p i are 
taken in independent simulations to explore the parame- 
ter space. The initial parameters adopted for our reference 
runs are summarized in Tab^ 

For planetesimal accretion to occur, a balance between 
different competing mechanisms is needed. For too small 
planetesimals, there is the risk that a too strong gaseous 
friction leads to a fast inward drift that prevents accre- 
tion of larger bodies far from the star. This is illustrated in 
Fig0J a population of 200m-sized bodies migrates towards 
the star in less than 2 x 10 3 years leaving the region be- 
yond 2 AU totally depleted of material. On the other hand, 
large planetesimals, that are less affected by gas drag, do 
not significantly migrate but their periastron alignment is 
weaker. If it gets too weak, then encounter velocities may 
not be reduced below the value that allows accretion into 
larger bodies. 
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Fig. 5. Eccentricity vs. semimajor axis plot after t = 
50000 years for lOkm-sized planetesimals. The dashed line 
is the corresponding average encounter velocity distribu- 
tion as a function of the distance to the star. Gas density 
is the same as in FigQ] 
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Fig. 6. Same as Fig0but for 50km-sized bodies 

In between these two accretion inhibiting cases, there is 
a planetesimal size-range where periastron alignment and 
eccentricity damping are efficient enough to reach very low 
Av without removing bodies from the system on a short 
timescale. As an example, for 10 km bodies the periastron 
alignment keeps the impact velocities below lOm/s within 
2.5 AU from the star (Fig(5J and the drift rate is also 
slow enough not to deplete this region of the disc within 
10 5 years. Fig|fj] shows the limiting case, i.e. the run with 
the biggest R p i (50 km) in which accretion is still possible 
at 2.1 AU. Beyond 2.2 AU the collision velocities begin to 
grow quickly as the secular perturbations are dephased. 

As the planetesimal accretion proceeds, larger bod- 
ies are formed and they still collide with the smaller 
ones that possibly still make up most of the mass of the 
swarm. It is relevant to verify whether the impact veloci- 
ties between large and small planetesimals still favour ac- 
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Fig. 7. Mutual relative velocities, at t — 30000 years, for 
a disc with 3 different populations of objects: 1 km, 10 km 
and 50 km in sizes. The solid line stands for Av among 
50 km bodies, the dashed line for Av between 50 km and 
10 km bodies and the dotted line for Av between 50 km 
and 1 km bodies. The gas disc has the same characteristics 
as in Figs l4l5l andlfil 
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Fig. 9. Impact velocities, at t — 15000years, on 500km 
bodies for 50km (continuous line), 10km (dased line) and 
1 km (dotted line) impactors. Same high density truncated 
disc as in FiglHl The absence of AV values beyond 2.5 AU 
for 1 km impactors is due to the fast removal of such small 
objects from these regions 
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Fig. 8. Same as FigEl(50km bodies at t — 50000 years) 
for a planetesimal disc initially truncated at 3AU. The 
gas density decreases radially as r~ 15 . 



cretion rather then fragmentation. The perihelion align- 
ment is indeed different depending on the size of the 
body. Small er bodies tend to align their perihelia to- 
wards 270° llMarzari and SchollL [2000) while larger plan- 
etesimals, less affected by the drag force, align to larger 
values. At the orbital crossing the impact velocities may 
thus be higher compared to those of equal size bodies. We 
performed additional simulations where we include dif- 
ferent size planetesimals. Fig0 shows the relative veloc- 
ities between populations of different size planetesimals. 
For 50 km target objects, colliding speeds significantly 
increase when the impactor sizes get smaller, especially 



in the otherwise low Av region below 2.5 AU. The Av 
might reach ~ 300m. s^ 1 when 10 km impactors are con- 
sidered, and even ~ 700m. s^ 1 for 1km objects. However, 
it can be shown that such impact velocities still lead to 
accretion. Adopting the collisional algorithm described in 
IPetit &: Farinellal l|l993j) , we find that a collision between a 
50 km and a 10 km object at 300 m.s -1 leads to net accre- 
tion of ~ 85 % of the impactor's mass, and that a collision 
between a 50 km and a 1km object at 700 m.s^ 1 leads to 
net accretion of ~ 95 % of the impactor's mass. 

Conditions slightly more favourable to planetary for- 
mation, i.e. allowing accretion of objects bigger than 
50 km, are met if we adopt a flatter gas density profile 
or if we truncate the disc closer to the star. In Fig|S]we 
show the relative velocity distribution in a disc of 50 km 
radius planetesimals cut at 3 AU from the central star and 
a gas density profile in r~ 15 granting a stronger drag force 
between 2 and 3 AU. Additional simulations show that un- 
der these conditions low Av are maintained for bodies up 
to ~ 200 km, while the inner limit of the low Av region 
never extends beyond 2.5 AU. 

For bodies in the 200-1000 km range, higher Av are 
obtained in the a < 2.5 AU region, especially for colli- 
sions with smaller planetesimals, which are probably the 
most frequent impactors on such bigger objects (Fig|!5J). 
However, such high Av impacts still lead to n et mass ac- 
cretion. Using again IPetit fc Farinellal l)l993|) algorithm 
for velocities obtained in Fig[5| on a 500 km target, one 
gets that 95% of the impactor's mass is reaccreted af- 
ter a 1200 m.s -1 collision with a 50 km impactor, 97% for 
a 1400 m.s -1 impact with a 10 km body and 98% for a 
1700 m.s -1 collision with a 1km impactor. So even con- 
sidering the large uncertainties on collision and reaccre- 



6 



Thebault et al.: The 7 Cephei System 



tion physics, it seems that we are far below the accretion 
threshold for 500km bodies in our "extreme" gas disc case, 
i.e. high gas density in the outer regions (flat radial distri- 
bution) and cut-off at 3 AU for the planetesimal swarm. 

Furthermore, even with lower gas densities, the final 
accretion phase of objects in the 200-1000 km range can 
be favoured by additional mechanisms that cannot be ac- 
counted for in our N-body simulations. Dynamical fric- 
tion produced by the gravitational interactions between 
planetesimals tends to prod uce equipartition of energy 
ijWetherill k. Steward Il989|) . As a consequence, larger 
bodies have lower random velocities favouring softer im- 
pacts. Self-gravity of planetesimals may also restore peri- 
astron alignment, as recently showed by Kokubo (personal 
communication). The scenario that comes out from our 
simulations can thus be summarized as following: Small 
to medium size planetesimals are strongly affected by gas 
drag and their relative velocities are kept low by the peri- 
astron alignment. When the gas drag weakens, additional 
dynamical mechanisms, such as dynamical friction and 
self-gravity, are still active and help large planetesimals 
to continue their growth. In case of massive protoplane- 
tary discs, gas drag can be able to do all the job, allowing 
the growth of planetesimals until they are big enough to 
even sustain high velocity impacts. 

It should also be noted that the high gas densities re- 
quired to maintain the periastron alignment up to large 
size planetesimals might not be unrealistic in the case of 
the 7 Cephei system. It is reasonable to expect that ac- 
cretion discs surrounding F2 stars, progenitor of K III 
giant stars as 7 Cephei are massive compared to the min- 
imum mass solar nebula. Even the density value for the 
gas adopted by l)Bodenheimer et ail [2000) for the disc 
surrounding 47 UMa, a GOV star, may be a lower limit 
for discs around F2 stars. 

An additional conundrum for planetary formation in 
7 Cephei is whether planetesimals could form in the outer 
regions of the disc perturbed by the gravity of a secondary 
star. How sensitive was the dust sticking mechanism to the 
gravitational pull of the star? It seems reasonable to as- 
sume that the sticking mechanism was more efficient in 
the inner disc, where the perturbations of the secondary 
star were weaker. The truncation of the planetesimal disc 
beyond a few AU from the star might thus be supported 
by the physics of planetesimal formation. At which radial 
distance did the truncation in the planetesimal population 
become significant? Only detailed models of planetesimal 
formation that include the gravitational effects of the com- 
panion star can answer this question. 

On the other hand, we should not forget that for too 
high gas densities, small bodies would spiral towards the 
star on a short timescale. The question is then: did ac- 
cretion proceed fast enough to prevent such a loss of 
small bodies where a large fraction of the mass lies? 
Planetesimal accretion in the 7 Cephei system thus re- 
quires a delicate balance between perihelia alignment at 
large sizes, fast spiralling for small bodies, and proper val- 
ues of the impact probability. 



A numerical model that includes all these effects at 
once is, at present, beyond computer capabilities. The 
main result of this section is that gaseous friction opens a 
window for starting planetesimal accretion within 2.5 AU 
from the star. 

3. Core formation by accretion of protoplanets 

In this section, we shall assume that planetesimal accre- 
tion could take place and lead to the formation of plane- 
tary embryos. This embr yo formation could have f ollowed 
the scenario described bv lKortenkamp et alJ l)200l[) . These 
authors have modeled planetesimal accretion in systems 
with a massive external perturber showing that the com- 
bination of gas drag, collisions, and secular perturbations 
of a massive external body favours orderly growth in the 
initial phases of planetary accretion followed by a phase 
of "type II" runaway growth. 

We here investigate whether giant impacts between 
such massive embryos in 7 Cephei can lead to the forma- 
tion of a planetary core whose orbit resembles that of the 
observed planet and whose mass is at least 10 M®, the 
mass required to trigger the final phase where the core 
accretes the re maining gas through rapid infall and forms 
a giant planet ( Po llack et We also evaluate how 

the timescale of the core formation depends on the initial 
mass of the embryos, on their number, and on their spatial 
distribution. 

We have simulated, within a full N-body model, the 
evolution of a population of protoplanets into a massive 
core by integrating the orbits of a swarm of planetary 
embryos distributed between 1 and 2.5 AU where, ac- 
cording to the previous computations, planetesimal accre- 
tion is mostly efficient. The orbital evolution of the em- 
bryos and their colli sions are comput ed with Chambers' 
Mercury code I Chambers et all l2002h which we modified 
to account for the large perturbations of the binary com- 
panion. The number and initial masses of the embryos 
are derived assuming a reference surface density of solid 
material at 2.1 AU ranging from 50 to lOOg cm -2 . As 
in the previous section, thes e values are derived from 
the iBodcnhei mer et 

ID (POOO) estimates for formation of 
the planetary companion of 47 UMa. The embryos initial 
masses range from a few Lunar mass to the mass of Mars, 
depending on the adopted value of the disc density. In 
all the simulations we assume that only a fraction of the 
mass of the solids in the disc has accumulated into plane- 
tary embryos while the remaining mass is still in smaller 
planetesimals. This percentage varies from 50% to 75% of 
the total mass in different models. We also adopted three 
different distributions of the protoplanets as a function of 
the distance from the star. This is motivated by the uncer- 
tainty about the details of the planetesimal accretion pro- 
cess. Runaway growth, type II runaway growth, orderly 
growth, or even oligarchic growth can occur depending on 
the delicate balance between the mas s and the velocity 
distribution. iKortenkamp et a (l200lh have shown that, 
for a particular binary system, type II runaway growth is 
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to be preferred. However, the type of growth strongly de- 
pends on the binary orbital and physical properties, and 
on the disc parameters like the density and mass distribu- 
tion. We prefer to cover a large spectrum of possibilities by 
assuming in our simulations different but equally possible 
initial conditions for the protoplanets: 

1) A population of embryos whose number and location 
is computed according to a superficial density a constant 
on average between 1 and 2.5 AU. 2) A a that decreases as 
\jr giving a constant mass for the embryos populating an- 
nular portions of the disc surrounding the star. 3) A fixed 
initial radial distance between the planetary embryos ex- 
pressed in Hill's radii. In some of the simulations we even 
include a " proto-core" , an embryo with an Earth mass 
located at 2.1 AU. Initially, all embryos have eccentrici- 
ties lower than 0.04 and inclinations lower than 1° with 
respect to the orbital plane of the binary system. 

Our simulations with an initial total embryo mass of 
25 M®, about 50% of the solid mass in a disc with a = 50 
g.cm -2 , fail to create a core of 10M® within lOMyr, the 
typical lifetime of a gaseous disc. The maximum mass 
of the core achieved in these simulations is 6 M® , which 
might not be enough to trigger the final gas accretion 
phase. If we increase the total embryo mass to 35 M®, 
a core of 8-10 M® can form in a few cases. All simulations 
with an initial mass ranging from 50 to 75 M®, compati- 
ble with a ~ 100 g cm" 2 , lead to the formation of a core 
with a mass up to 20 M® within 10 Myr. However, our 
simulations all show that the core always ends up within 
1.5 AU from the primary star, while the observed planet is 
at about 2.1 AU. Even when including in the initial proto- 
planet population a bigger "proto-core" of 1 Earth mass 
at 2.1 AU, the final accreted core always ends up between 
1 and 1.5 AU. Two distinct mechanisms might account for 
this outcome: 1) In some cases a core might begin to form 
in the outer regions but it migrates inwards due to plan- 
etesimal scattering reinforced by the gravitational forces 
exerted by the binary star which excites the embryos' ec- 
centricities fFig llOfl . 2) In some other cases a core begins 
to form at around 1 AU. It grows at a faster rate because 
of the shorter Keplerian orbital period of bodies in in- 
ner orbits. This core is able to accrete large protoplanets 
in outer orbits, since the orbital eccentricity is large due 
to the companion's perturbations, and the final result is 
a core around 1.5 AU and almost no remaining material 
beyond this position. This result hold even when a large 
"proto-core" is initially placed in the outer disc fFig lll(l . 
These two mechanisms are peculiar of protoplanetary ac- 
cretion in binary star systems, where large eccentricities 
among the embryos are excited by the gravitational pull of 
the secondary star. In a few simulations the 'proto-core' 
is ejected out of the system before it reaches 10 M® , since 
it is close to the border of the stability region. 

In Fig 1 121 we show the outcome of a simulation with 
an intermediate value for the total embryo mass (75 M®) 
and a proto-core initially located at 2.3 AU. The proto- 
core grows faster than nearby protoplanets but it migrates 
inward due to the scattering of the other bodies. It settles 
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Fig. 10. Migration of the growing core due to scattering 
of protoplanets. The initial disc extends from 1.5 to 2.5 
AU and has a total mass of 75 M® . 
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Fig. 11. Evolution, in the (M,a) plane, of a 1 M® proto- 
core initially put in the outer disc and of the largest em- 
bryo growing in the inner region. The rapid growth of the 
inner embryo results in the capture of the outer proto-core 
(see text for details). 

at about 1.4 AU and its mass reaches almost 10 M®. The 
final giant planet is then expected to orbit closer to the 
primary star than the observed planet in 7 Cephei 

4. Discussion and conclusions 

The previous results show that planetesimal accretion in 
7 Cephei depends on a delicate balancing between gas 
drag and secular perturbations by the secondary star. 
If this balancing is met, then type II runaway growth 
l|Kortenkamp et all l200lf) might possibly lead to the for- 
mation of massive planetary embryos within 2.5 AU from 
the central star. In the subsequent phase when giant im- 
pacts between the embryos build up a massive core, the 
major problem is not the timescale but the final location 
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instability scenario but this scenario remains 

to be quantitatively tested for the 7 Cephei system. 
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Fig. 12. Snapshots of protoplanet accretion in the 
7 Cephei system. The initial mass in the embryos is 75 
M s and the superficial density a decreases with the ra- 
dial distance as 1/r. 



of the planet, which is always well inside the actual po- 
sition of the observed planet. Several explanations might 
be proposed to account for this discrepancy. A first hy- 
pothesis involves the evolution of the whole system. It is 
possible that the distance between the two stars was larger 
when the planet formed and that after the planet forma- 
tion additional mechanisms pushed the secondary star on 
an inner orbit. This could be the case if the 7 Cephei sys- 
tem was born in a clustered environment, where close en- 
counters with other young stars may caus e perturbation s 
of the binary orbit that tend to shrink it jHeggiel Il985[) . 
A different way to reduce the orbit of a binary system is 
related to the possibility that, originally, the system was 
triple or more. The ejection of one or more stars causes 
a transfer of binding ene rgy and an eve ntual reduction 
of the binary separation (jReipurtbl l2000h . A more com- 
plex mechanism is related to the formation of more than 
one giant planet around 7 Cephei . If the circumstellar disc 
around the main star was significantly more massive, with 
a superficial density of solids higher than 100 g cm -2 , it is 
possible that two or more giant planets formed around the 
main star. Mutual scattering among these planets ejected 
one, two or more of them out of the system leaving a 
single planet in the observed orbit l|Marzari et all l2004|) . 
Another and more radical solution would be to renounce 
the core-accretion model in favour of the alternative disc 
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